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Submillimeter Continuum Properties of Cold Dust 
in the Inner Disk and Outflows of M 82 

Lerothodi L. Leeuw^'^'^ and E. Ian Robson^ 
ABSTRACT 

Deep submillimeter (submm) continuum imaging observations of the starburst 
galaxy M82 are presented at 350, 450, 750 and 850 /xm wavelengths, that were 
undertaken with the Submillimetre Common-User Bolometer Array (SCUBA) 
on the James Clerk Maxwell Telescope in Hawaii. The presented maps include a 
co-addition of submm data mined from the SCUBA Data Archive. The co-added 
data produce the deepest submm continuum maps yet of M82, in which low- 
level 850 /im continuum has been detected out to 1.5 kpc, at least 10% farther in 
radius than any previously published submm detections of this galaxy. The over- 
all submm morphology and spatial spectral energy distribution of M 82 have a 
general north-south asymmetry consistent with Hq, and X-ray winds, supporting 
the association of the extended continuum with outflows of dust grains from the 
disk into the halo. The new data raise interesting points about the origin and 
structure of the submm emission in the inner disk of M 82. In particular, SCUBA 
short wavelength evidence of submm continuum peaks that are asymmetrically 
distributed along the galactic disk suggests the inner-disk emission is re-radiation 
from dust concentrations along a bar (or perhaps a spiral) rather than edges of 
a dust torus, as is commonly assumed. Higher resolution submm interferome- 
tery data from the Smithsonian Submillimeter Array and later Atacama Large 
Millimeter Array should spatially resolve and further constrain the reported dust 
emission structures in M82. 
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Introduction 



Massive ejections of gas and dust originating in galactic nuclei have been observed at 
scales of a few kpc in optical emission lines, submillimeter (submm) molecular lines, mid- 



(e.g., 


Watson et a 


. 1984 


Eneelbracht et al. 


2006 





198J; iHeckman et all Il990l : IPevine fc Ballvl Il999l : iHoopes et all 12005 



One explanation for the outflows is that a high supernova rate in 
the galactic nucleus heats up the surrounding gas to high temperatures with sound speeds 
exceedi ng the escape velocity of the galaxy, creating a wind that expands outward from the 
galaxy (jChevalier &: Cleggi Il985l ). The wind entrains cosmic rays, warm and cool gas, as 
well as cool dust, making the outflow directly visible in many wavebands. The outflows are 
usually oriented along the minor axes of the galaxies and are thus most easily observed in 
edge-on galaxies. In local starbursts and high-z Lyman Break galaxies with high enough 
global star-formation rate per unit area, the superwinds are common and responsible for 
expelling metals from these galaxies and enriching the inter-galactic medium (I GM) and 
therefore have implicatio ns in the evolution of galaxies and the IGM (see, e.g, 
2003: Veilleux et al. 2005. for recent reviews). 



Heckman 



M82 (NGC3034) is a nearby and popular object in which to investigate the physical 
association between galactic nuclei and large-scale outflows. The galaxy is edge-on with 
an inclination of about 10° at position angle 72° and is classified as Irr ll. At an estimated 
distance of 3.63 Mpc (as determined for M81 by iFreedman et al.l (119941 )). M82 has optical 
dimensions of 11'. 2 x 4'. 3, i.e. ~ 11.8 x 4.5 kpc (image scale is ~17.6 arcsec"^). The nuclear 
region, within 4' x 2' about the major axis of the galactic disk, has numerous point sources or 
emission concentrations, some of which originate from supernovae and massive star clusters 
and have been detected from the X-ray to radio wavebands. Layers of dust filaments laden 
these inner regions producing severe optical extinction and copious infrared to submm re- 
radiated emission. 

New Hubble Heritage Team optical images obtained with a deep 6-point mosaic in B 
(0.45 /im), V (0.55 /xm), I (0.81 /im), and H^ (0.65 /xm) filters of the Advance Camera for 
Surveys (ACS) on bo ard the Hubble Space Telescope (HST) exhibit detailed, filamentary 
outflows of the M82 (IMutchler et al.l 120071 ). especially in H^. As described above, it is 
thought that this outflow is being driven by the copious formation of massive stars (or a 
starburst) and subsequent explosions of supernovae. The starburst outflow not only provides 
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the ejection mechanism for the material from the galactic nucleus, but also heats the gas 
and ionizes the hydrogen, causing it to glow with the red light of emission line. 

These new optical images as well as earlier detailed natural-color comp osite images of 



M 82 obtained with HST (see, e.g.. lde Grijsll200ll ) and the Subaru Telescope (lOhyama et al. 



2OO2I ) show more than 100 compact groupings of about 10^ stars in very bright star clusters 
sprinkled throughout the galaxy's central region, prominent dust lanes that crisscross the 
disk, knotty filaments of ionized gas that have rich nebular spectra that are not especially 
enriched in nitrogen, hydrogen gas in a strong galactic wind that is clearly below the galactic 
center and to the right of the central region, along with many other regions of varying star- 
formation environments in the nuclear parts of this galaxy. The huge clusters of massive 
stars, numerous X-ray and radio detected supernovae, gas concentrations, optically-dramatic 
dust filaments, galactic winds and other active nuclear features have been attributed to a 
large burst of star formation 10'' to 10^ years ago, that was probably triggered by a tidal 
inte raction with the nearby spiral galaxy M81 and dwarf starburs t galaxy NGC3077' (see , 
e.g.. lYun et al.l (119941 ) for evidence of H I tails linking the three and lForster Schreiberl (120001 ) 
for a recent review). 

The proximity of this galaxy makes it possible to observe the region of interaction be- 
tween the star-formation regions and the halo, including expanding shells or bubbles and 
"chimneys" that are producing a clearer picture of the l ocalized driving mechanis ms for 



the outflows (e.g., iHeckman et al.lll990l : IWills et al.lll999l : IWestmoquette et all 120071 ). The 



proximity also allows the detections of low-level emission in the halo and consequently the 
deter mination of the amounts and possible origin of material that results in this emission 



[e.g., iSeaquist fc ClarkI I2OOII : lEngelbracht et al.l 120061 ). Studying the contents and interac- 
tions between the star- formation regions and the halo is important for understanding their 
role in the evolution of M82 and may provide clues to general galaxy evolution as well as 
details of the composition of the intergalactic material. 

This paper focuses on the Submillimetre Common-User Bolometer Array (SCUBA) 
maps of the copious submm re-radiated emission that results from the dust-laden, star- 
forming disk, as well as large-scale, low-level emission that is associated with the outflows 
in the halo of M82. Submm continuum observations of the dusty central regions in M82 
were previously obtained with the submm co n tinuum rece iver UKT 14 on the James Clerk 



Maxwell Teles cope (JCMT) bv iHughes et al.l ( 



1990 



199J) and later with SCUBA, also on 



the JCMT, bv lleeuw etall Jl999h and klton et airjl999h . The current study was intended 



to extend these previous imaging submm observations in spatial extent, sensitivity and to 
all the four submm wavelengths available with the SCUBA array (see Section [2]). The study 
includes a co-addition of data mined from the SCUBA Data Archive. 
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The co-added data produce the deepest submm continuum maps yet of M82, in which 
low-level emission is detected out to 1.5 kpc for the first time in the submm continuum of 
this galaxy. The deep maps are used in a detailed morphological study of the nuclear and 
large-scale detections (see Sections [3] and ED, including a focused comparative analysis to op- 
tical (see Section [3l2!) and high-resolution CO (1-0) (see Section ISTT!) morphology. The maps 
are also used in the computation the first submm spatial spectral energy distribution (SED) 
of separate locations in the nuclear star- forming region of M82 (see Section [5]). These obser- 
vational results are used in the discussion of the origin and structure of submm continuum 
morphology and spatial SED of M 82 and reviewed in the context of relevant interpretations 
by other researchers, including those who use data from other wavelengths (see Section [5]). 
In particular, (a) the commonly assumed interpretation that the double emission peaks that 
are seen in the mm to infrared continuum are due to emission from the edges of an inclined, 
dusty-molecular torus is challenged (see Section EH]), (b) an analytical review of CO results 
is undertaken to assess if CO emission may significantly contaminate the continuum observed 
in SCUBA filters (see Section Hj), and (c) a morphological comparison is conducted to check 
whether the loc alized outflows tha t are reported in hig h res olution radio, CO, an d SiO m aps 
respectively bv Iwills et aP Jl999h . IWeiB et all JlQQoh . and iGarcia-Burillo et all J200lh can 
be seen in the SCUBA maps (see Section 1531) . Finally, the overall implication of the results 
are discussed and possible future work is outlined (see Section [6]). 



2. Observations 

SCUBA 850, 750, 450 and 350 /xm imaging observations of M82 were obtained with 
the telescope poin ted at the 2.2 /xm infrared nuclear peak of the galaxy using positions from 



Dietz et al.l (119861 ). Jiggle mapping observations were conducted with the secondary chopping 
in azimuth at 7.8 Hz and with a throw of 120". The imaging observations employed the 
common 64-point jiggle pattern with a 3 " offset between each position, giving fully sampled 
images with both arrays. 

Because M 82 has been a popularly observed source with SCUBA, additional maps that 
were obtained with a chop throw of 120 " by other observers were mined from the SCUBA 
Data Archive in order to co-add the related data and maximize the signal-to-noise in the final 
maps. Data sets from the SCUBA Archive were separately flux calibrated and corrected for 
JCMT pointing errors and then co-added, with each observation being weighted according 
to its relative integration time and the noise in the map. The 450-850 /xm dual mapping 
wavebands have been used more commonly than the 350-750 /xm combination, and therefore 
the archival maps constitute about 85% and 35% of the respective co-added, total-integration 
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time for the 450-850 fim and 350-750 fim dual maps. 



The imagin g data analysis was u ndertaken using the dedicated SCUBA data reduction 
software SURF flJenness et aPllQQsh . as well as KAPPA, GAIA and CONVERT software 
packages provided by the Starlink Projeclll]. The data reduction consisted of first fiatfielding 
the array images, and then correcting for atmospheric extinction. Next, pixels significantly 
noisier than the mean were blanked-out; and, after initial inspection of raw images, pixels 
containing relatively little flux from the source were used to correct for correlated sky noise 
in each individual jiggle-map. 

The atmospheric opacity, r, was determined from skydips made with SCUBA at intervals 
during the observations. Were SCUBA-skydip measurements were not available, the r at 
the SCUBA filters was extrapolated from the continuously measured r at 225 GHz, obtained 
courtesy of the Caltech Submillimeter Observatory (CSO) radiometer and using relations 
listed in Table [H At the JCMT, these relations are empirically derived and periodically 
updated and impro ved as more data, e specially since the commissioning of SCUBA, have 
been obtained (e.g., Archibald et al.| [2002lR. All the data were calibrated using instrumental 
gains determined primarily from nightly beam maps of Mars and Uranus, or, alternatively, 
the JCMT secondary calibrators. 



Table 1: The CSO relations used on the data presented in this paper 



T850 


= 4.3 X (rcso - 


- 0.007) 


''"450 


= 23.9 X (rcso 


-0.01) 


''"450 


= 6.5 X (rgso - 


0.03) 



3. The General Submm Continuum Morphology 

Figure [1] shows the 850, 750, 450 and 350 fim co-added maps of M82 that were obtained 
in jiggle mapping observing mode, respectively at 14". 5, 11". 4, 8". 5, and 6". 7 resolution and 
about 8, 100, 225, and 650mJy/beam sensitivity. The 850 and 750 /im images have a single 
emission peak that is centered about 9" west of the galactic nucleus, while the 450 and 
350 /im maps have two emission peaks centered about 10" and 6" respectively east and west 



^The Starlink Project is run by the Council for the Central Laboratory of the Research Councils on behalf 
of the Particle Physics and Astronomy Research Council of the United Kingdom. 

^This latest documentation is available at the Joint Astronomy Centre web site, 
http: / / www.jach.hawaii.edu| 
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of the nucleus along the galactic disk. The peak in the west is slightly elongated along the 
galactic disk and is brighter than the eastern peak, showing east-west asymmetry about the 
nucleus. In the 750/Ltm image, the single emission peak seen in the 850|Ltm image begins to be 
resolved-out into the double peaks seen in the 450 and 350/im images (e.g. Figures [U bottom 
panel). This is expected, as the 750/im observations have a resolution that is intermediate 
between that of the 850 and 450|Um. When the maps at shorter wavelengths (e.g., 350/um) 
are smoothed to resolutions similar to those at the longer wavelengths (e.g., 850/im), the two 
submm peaks that are resolved at the shorter wavelengths become visible just as one peak, 
showing that indeed the different high-brightness morphologies depicted in the maps are 
due to the respective resolutions at different wavelengths. The overall extended morphology 
has an elliptical shape with the major axis position angle of 72°, i.e. roughly the same 
as the galactic disk in the nuclear reg ion. This general mor p hology is sim i lar to previous 
continuum observat i ons at mm (e.g.. 



flHughes et all Il994l : iLeeuw et all Il999 



Kuno fc Matsuol Il997l : iThuma et aP boooh . submm 



Alton et al. 19991 ) and mid-infrared ( Telesco et al. 



199lD wavelengths, as well as to CO line transitions (e.g., iNakai et al.l 119871 : iThuma et al. 
2000, ) observations of comparable resolution. 



It is noted that although the maps p resented in Figure [J have similar featu r es to those 
in published submm continuum maps by iHughes et al. Jl994h and lAlton et~aD Jl999h . the 
very sensitive 850 /im map presented here also shows that the 850 ^m emission (that is 
expected to be from cold dust) extends by at least 10" (~ 176 pc) radius farther into the 
halo than detected by those authors. Furthermore, Figure [2] depicts the 850 /im continuum 
emission and integrated C0(2-l) li ne intensity maps re spectively presented in Figure [T] of 
this paper and Figure 5 of work by iThuma et al.l (120001 ) to sh ow that the 850 emission 
is as extended as the C0(2-l) emission, contrary to claims by lThuma et al.l (120001 ) that the 
CO (2-1) emission is more extended than the cold dust emission of M82, and therefore made 
the galaxy exceptional in this regard. 



3.1. Submm vs. High-Resolution CO (1-0) Morphology 

Figure [3] shows black contours of the integrated-intensity-CO (1-0) data, th at were ob- 
tained with the Berkeley Illinois Maryland Association (BIMA) interferometer by lShen fc Lo 
(119951 ). overlaid on the SCUBA 450 and 850 fim continuum images by aligning the sky co- 
ordinates in the SCUBA maps to those of the BIMA ones. The BIMA maps are at resolution 
2". 4 X 2". 6 and are plotted to investigate the spatial correspondence between the submm 
continuum morphology and the high-resolution CO features. The top panels in Figure [3] are 
plotted at full resolution of the BIMA CO contours and the bottom panels are with the CO 
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contours smoothed to beam sizes similar to the SCUBA-850 and -450 fim beams (~ 14". 5 
and 8". 5, respectively). Because dust often occurs mixed-in with gas in star-forming regions, 
the BIMA maps are expected to give an indication of how dust emission may appear at 
higher resolution and perhaps also some insight about the structure seen in the SCUBA 
images. The BIMA maps probably provide the best observational, high-resolution evidence 
of structure in the inner disk that is associated with dust, as submm continuum observations 
of dust in M82 are not currently publicly available at a resolution higher than 6". 7, that is 
obtainable with SCUBA at 350 fim. 

The first obvious difference between the SCUBA and high-resolution-BIMA maps is 
that the dust emission peak indicated by a red contour is resolved into two peaks centered 
about 9" apart in the CO maps. Of these two CO western peaks, the one closer to the 
nucleus is co-spatial with the 450 fim (red contour) and slightly east of the 850 fim (red 
contour) western peaks. The CO peak farthest from the nucleus is actually the brightest 
of all the CO peaks and lies on the eastern edges of the submm 850 and 450 /im high- 
brightness lobes, not co-incident with the brightest submm peaks. Distinct from submm 
continuum, the lower level CO emission near the very western CO peak fans-out westerly 
at position angle 85°, diverting from the 72°-position-angle of the submm continuum lobes 
and disk along the major axis, as well as that of the CO within a 10" radius from the 
galactic nucleus. These differences in the CO and dust features are evident in both the full 
resolution and smoothed CO maps plotted in Figure [3], suggesting that although the dust and 
CO generally appear mixed in this star-forming region, in fact there are differences in their 
spatial distributions. The different concentrations are most probably locations of varying 
gas-to-dust densities or star-formation environments. Consistent with the finding here o f 



different CO-dust con centrations, evidence of a star-formation history (e.g.. Ide Grijsll200ll ) 



and gas density (e.g., iPetitpas fc WilsonI |2000( ) that clearly varies from the east to west of 



the galaxy has been reported in M82. 

Higher resolution submm observations of dust in M82 that should be possible with 
the Smithsonian Submillimeter Array (SMA) and later Atacama Large Millimeter Array 
(ALMA) will provide direct observational evidence to further test how different the CO and 
dust emission trace each other at the small scales shown in the BIMA maps. Those future 
observations will also test if dust emission has more complex morphology than has currently 
been detected, as is suggested by the increasing structure that is seen in the SCUBA maps 
going from low to high resolutions. Further discussion of the origin and structure of the 
submm emission peaks is detailed in Section O 
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3.2. Submm vs. Infrared and Optical Morphology 



The spatial investigation of submm versus infrared and optical morphology in M82 
is important because complex optical morphology that is s een in this galaxy, with vis ual 



extinction values [Ay) that range from about 3 to 25 (e.g., lAlonso-Herrero et al.ll200ll ). is 
thought to result from obscuration of optical light by large, cold dust grains that are heated 
by stars (among other things) a nd re-radiate in the infrared to subr nm wavelengths. Evidenc e 
of star-formation history (e.g.. Ide Griisll200ll ). gas density (e.g.. iPetitpas fc WilsonI l2000l ) . 
and submm and CO emission peaks (this work, e.g. Section 13.11) that clearly varies from 
the east to west of the galaxy strongly suggest that any associated dust lanes must vary not 
only in their geometric structure but also in their heating mechanisms and composition. 



3.2.1. Strong Optical- Oh suration Patches that Correspond with Submm and CO Peaks 



T he panels in Figure HI show the S-band maps of M 82 obtained by the Hubble Legacy 
Team ( Mutchler et al.ll2007l ). overlayed with the 450-)um continuum -emission contours shown 
in Figures [T] and CO (1-0) interferometry data by lShen fc Lol (119951 ). The overlays were made 
by aligning the sky co-ordinates of the Hubble maps with those of the SCUBA and BIMA 
ones. The respective resolutions in the Hubble, SCUBA, and BIMA maps are ~ 0.09" 
~ 8". 5, 2". 5. For best contrast, the optical intensities in this figure are inversely plotted, 
and therefore the light patches are strong extinction features. The panels show that all the 
submm and CO peaks of M 82 are spatially coincident with very strong optical-obscrution 
patches seen with in the intense sta r-forming region about 30" x 15" of the galactic nucleus 
and 2.2 /zm peak (jPietz et al.lll986l . marked by a cross in the figures). The brightest submm 
peak that is about 9" southwest of the nucleus and the CO peaks in this region stretch along 
the same direction as prominent optical dust lanes that lie east-west along the galactic disk 
major axis and co-spatial to these dust lanes out to about 15" from the galactic nucleus. 
Interferometery maps of the central CO(l-O) peaks shown in the left panel of Figure H] also 
demonstrate the high-resolution, spatial co-incidence between the CO (1-0) peaks and very 
dense optical-obscrution patches about the nucleus. The most western and brightest CO 
peak also co-coincides with a dense optical-obscuration patch. 

Within a radi us of about 1 0" ab out the galactic nucleus, or about two arcseconds north 
of the 2.2 fim peak lDietz et al.l (119861 ). there is intense 5— band optical emission. This region 
is between the two submm peaks and thus has relatively low submm intensity or, if indeed 
the submm emission is from cold dust re-radiation, low density or heatin g of cold dust. This 
regio n also has mid-IR emission indicative of star-formation clusters (ILipscy fc Plavchan 
20041 ). though of less mid-IR brightness and lower mid-IR color temperature than the star- 
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formation clusters southwest of this position, i.e. at the location of the southwestern submm 
peak. The remarkable spatial co-incidence between the submm as well as high-resolution CO 
peaks and very dense optical-obscuration patches about the nucleus, and the co-incidence 
of intense S— band optical emission with location of relatively low submm emission or cold 
dust column density, support suggestions in this paper that the submm emission in the inner 
disk of M 82 originates from re-radiation of dust concentrations or clouds of physically dif- 



ferent star-for mation environments (e.g., lAchtermann fc Lacylll995l : iForster Schreiber 



2000 



de Grijsll200ll ). rather than the commonly assumed interpretation of a dusty torus about the 



nucleus (see Section ISTTl below) . 



In color images of M82 (see, e.g., Ide Grijsll200ll : IWestmoquette et al.l 120071 ). the opti- 
cal emission associated with submm peaks and high-brightness, diffuse re-radiation in the 
inner disk of this galaxy has a blue-brownish hue, strongly suggesting that hot, young, 
blue stars are the main heating source for the dust in this inner region. The young stars 
also produce infrared emission and emission lines associated with intense star-formation this 



galaxy (e.g.. lAchtermann fc Lacvlll995l : lAlonso-Herrero et al.ll200ll : iLipscy fc Plavchanll2004 



Westmoquette et al.l 120071 ). 



3.2.2. Strong Optical- Oh suration Patches with No Corresponding Submm and CO Peaks 

Although all the submm and CO peaks correspond to optical features in the inner disk 
of M82, as described above, the contrary is not true; i.e, many prominent dust, optical 
filaments, clouds, and lanes have no clear corresponding submm emission counterparts. In 
particular, the submm emission is basically smooth at the locations of 1) north-south fila- 
ments that run below, through, and flare above the southwestern submm peak that is about 
8" from the galactic nucleus; 2) east-west dark lanes that run west along the galactic major 
axis and continue through to about 35" west of the galactic nucleas (or 25" west of the 
southwestern submm peak); 3) huge, dramatic complex of optically-obscuring clouds, fila- 
ments and lanes that extends the entire minor axis of M82's disk and covers an area greater 
than a diameter of 25" south of the disk just east of the eastern submm peak that is about 
10" from the galactic nucleus; and 4) light, optically-obscu ring clouds an d filaments in a 



'low'-extinction region known as a starburst-remnant (e.g.. Ide Grijd l200ll ) that lies about 
30" to 60" northeast of the galactic nucleus. In another region of low-optical extinction and 
radiation that is about 60" to 120" northeast of the galactic nucleus, submm emission has 
currently not been detected where very light optically-obscuring filaments are evident. 



The lack of correspondence between the very dark optical clouds, filaments, and lanes 
with bright submm emission suggests that these optical features are due to obscuration by 
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cool dust grains that are on the near side of the galaxy and at large distances from the nuclear 
region, where they are heated by a very dilute stellar radiation field. These foreground dust 
clouds evidently have enough column densities to obscure optical light in the line of sight; 
however, they only re-emit very low-level emission that shows no striking features in the 
current submm maps. That the complex optical morphology seen in Figure H] is primarily 
due to foreground dust is supported by the fact that the obscuration is more dramatic in 
the shorter wavelength 5— band than longer /—band images (L. L. Leeuw et al. 2009 in 
preparation). Obscuring dust is expected to be optically thicker at the shorter wavelength 
and therefore cause more optical extinction and thus appear more prominently at the shorter 
wavebands. 

One explanation for the submm low-level continuum having a relatively smooth mor- 
phology is that, because this emission is optically thin, the detected radiation at a particular 
submm wavelength represents the total emission from the entire galactic column of dust in 
the line of sight. This is different from the optical morphology of dust because the dust is 
typically seen obscuring stellar light, and therefore only the dust in certain spatial strati- 
fications (usually the foreground) of the line-of-sight columns is observed. In other words, 
morphology due to spatial-depth or stratification of similar dust grains that are heated by a 
dilute radiation field is most often more obviously seen in optical obscuration than in submm 
emission. These morphological effects will of course depend on the sensitivities and resolu- 
tions of the instruments used. In a low optical extinction and radiation region about 60" 
to 120" northeast of the galactic nucleus, for example, the lack of any submm detection to 
date may simply because current continuum instruments have not been sensitive enough to 
easily detect low-level, dust re-radition that my correspondence to low-level extinction and 
stellar heating in those regions. The relatively larger SCUBA beam could also beam-smear 
and thus erase small features detected in the higher resolution optical maps. 



3.2.3. The Outflowing Wind 



An alternative explaination for the submm morphology of M82 being smooth as op- 
posed to distrupted like the optical morphology depicted in Figure H] is that the submm 
low-level emission is primarily due to dust entrained in outflowing gas and physically dif- 
ferent from the optical extinction features that don't have any currently detected submm 
counterparts. Figure E] shows the low- and high-br ightness features of the Hq, maps obtained 
by the Hubble Legacy Team (IMutchler et al.ll2007l ). overlayed with SCUBA 850 and 450-/um 
continuum-emission contours as shown in Figures [H For best contrast, the optical intensities 
are inversely plotted and the light patches (e.g. across the center of the image in the right 
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panel) are foreground extinction features. The emission is plotted saturated to high-light 
the large-scale and base of the outflowing wind, so not all the obscuration patches that 
criss-cross optical maps of M 82 are depicted here. The overlay in the left panel demonstrates 
the spatial co-incidence between the large-scale low-brightness Hq, emission and 850-/xm con- 
tinuum north and south of the disk, while the right panel depicts the origin of the large-scale 
Hq emission in the intense star-formation inner disk of M82 about the submm emission 
peaks (especially near the s outhwestern peak ), where a ~ 130 pc expanding "su perbubble" 
has been discovered in CO (IWeifi et al.lll999l ) and ionized gas (IWills et al.lll999l ). 



Recent large-scale, high- resolution OVRO observations by lWalter et al.l (120021 ) detected, 
resolved molecular CO(l-O) streamers in and below M82's disk that have different kinemat- 
ical signatures to its outflowing gas. Some of the streamers are well correlated w ith optical 



obscuration features and form the basis of some prominent tidal HI features (lYun et al. 



19931 ) that are thought to provide evidence that the gas within the optical disk of M82 is 



disrupted by the interaction of M 82 with M 81 and likely triggers the starburst activity in 
M82's center (jWalter et al.ll2002l ). The detection of resolved mid-IR to submm emission that 
corresponds to the optical and gas streamers and is perhaps a physically separate component 
to the outflowing gas and dust in M 82 should be possible with sensitive and high-resolution 
mid-IR to submm ima ging instruments u s ing Sp itzer, ALMA and the SMA. New Spitzer ob- 
servations reported by lEngelbracht et al.l (120061 ) did indeed detect extended mid-IR emission 
not only in the onflow of M 82 but also in its halo. The extended halo mid-IR emission could 
be from ma terial ejected into the ha lo by the outflowing wind or from the interaction of M 82 
with M81 (lEngelbracht et al.ll2006l ) . Future observations with these sensitive instruments 
and their detailed data analysis have the potential to 1) directly uncover the disruption 
of cold (and warm) dust distribution by the interaction of M82 with M81, 2) decompose 
submm dust emission in the disk, outflows, and streamers of M 82 and better constrain the 
properties of cold dust in these seperate components (see Section ESD, and 3) elucidate the 



role or consequence of the dust in the interaction of M81 and M82 (e.g., lYun et al. 
and any connected triggering and evolution of the star-forming in M82 (cf. 



1993h 



Walter et al. 



20021 ) . and the re-processing of galactic dust in general. 



4. Possible CO Contamination of the Submm Continuum? 

It is worth noting that the dust morphology that is mapped in the submm continuum 
from M82 may be substantially enhanced by CO emission from this galaxy. A recent flux 
comparison between CO (3-2) emission and 850 /xm continuum in M82 showed that CO 
makes a 47% (i.e. high) contribution to the integrated continuum in this SCUBA band 
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( ISeaquist fc Clarkll200ll ). F rom analyzing collate d CO (4-3) observations together with those 
of CO at lower transitions, iGuesten et al.l (119931 ) concluded that the line strengths in M 82 
increased as one went to higher transitions, indicating that the higher transitions must 
provide significant cooling in the galaxy. All SCUBA bands have roughly the same widths, 
i.e. 30 GHz, and therefore the CO contribution to the higher frequency continuum would 
be expected to b e equa lly or more significant than that reported for the 850 /im band by 
Seaquist fc ClarkI (120011 ). However, because the submm continuum in M82 has a thermal 
spectrum (see Section 15.21) . and therefore the submm fiuxes increase with frequency, the 
CO contribution to the higher frequency continuum may be less than the estimates for 
the 850 fim band. For SED analysis in this paper (e.g.. Section [5.21) . the CO percentage 
contribution to the measured fiux is assumed to be the same across the SCUBA bands and 
no correction for it is made in the presented data. 

Observations using a new high frequency, Fabry- Perot spectrometer on the JCMT have 



lead to clear detections of the high transition CO(7-6) in M82 and NGC253 ([Bradford et al. 



1999 ). the first such detection in any extragalactic sources. The analysis by iGuesten et al. 



(I1993I ) and the detections of CO (7-6), whose transition line lies in the 450 /im filter bandpass, 
suggest that other higher transition lines, such as the 13CO(8-7) line that lies at the centre 
of the SCUBA 350 fim band, could be very strong in M82, supporting the proposition above 
that the high frequency SCUBA images may have significant contribution from CO. In this 
light, the morphology seen in the SCUBA images is a direct probe of the galactic cooling 
and the general interactions of active star formation and the ISM in M82. 

Although it is not obvious if the CO contribu tion to the higher freque ncies of SCUBA 
will be less than or as significant as estimates by ISeaquist fc ClarkI (120011 ). it is clear that 
the CO contamination to higher-frequency continuum warrants investigation. Future work 
on this galaxy will attempt to acquire the data of the CO lines in the 450 and 350 /im 
bands and make a quantitative comparison of these data in order to determine the possible 
contributions of CO to the high-frequency-SCUBA data. Such work is important (among 
other things) in the determination and interpretation of submm SED and thus the nature of 
dust emission in M82 (see Section \52\ below) . 



5. Origin of Submm Continuum and Spectral Energy Distribution 



5.1. Source and Structure of the Submm Continuum in the Inner Disk 



The radiation frorn M 82 at the radio (e.g.. lSeaauist &: Odegard 



mm (e.g., 



Hughes et al.lll990l ) and submm-to-infrared (e.g. 



199 



Telesco et al 



Wi 



1991 



Is et al. 



Hughes et al. 
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1994 : lAlton et al.lll999l ) wavelengths is respectively dominated by synchrotron emission from 
supernovae, free-free emission from ionized gas, and thermal re-radiation from dust heated 
by young stars. In this light, the double peaks seen in the mm to infrared continuum have 
commonly been interpreted as due to emission from the edges of an inclined, dusty-molecular 
torus, that - as a result of their geometry on the plane of the sky and op tically-thin nature o f 
radiation - have relatively high optical-depths in the line of sight (e.g.. iHughes et al.lll994l ). 
In Figure [H and other mm-to-infrared maps of similar or worse resolution (such as those 
from IRAS), the double peaks are not resolved and appear as a single, elongated lobe that 
is brightest in the southwest. Like the galactic disk, the lobe (or peaks when resolyed-ou t) 
has a position angle of roughly 72° and - in the tori interpretation (e.g., Shen fc Lo 19951 ) - 
an inclination of ~ 10°. 

The peaks of emission seen in the mm to infrared have alternatively been interpreted sim- 
ply a s dust and molecular c oncentrations along the ga lactic disk, perhaps in a bar structure 



le.g- 



Neininger et al.lll998l : IWestmoquette et al.ll2007l ) that m ay have an expanding "super 



bubb le" of gas centered at supernova remnant 41.9+58 (e.g., IWeifi et al.lll999l : IWills et al. 



19991 ). This interpretation is supported by at least three reasons. First, the east- west, dou- 
ble peaks have now been seen in maps of both optically thin and thick CO emission (e.g., 
Petitpas &: Wilsonll2000l ). This is reasonable if the emission is from a structure that consti- 
tutes concentrations or clouds of dust but is contrary to what is expected if the emission 
is from a stru cture with tori ge ometry that, like the galactic disk, is thought to be highly 



inclined (e.g., IShen &: Lol 119951 ). For optically thin radiation, it will be possible to detect 



emission from the inner parts of the imaged structure, and either dust clouds or indeed edges 
of an inclined torus would manifest as regions of relatively higher optical-depth or brighter 
opt ically-thin emiss i on in the li ne of sight (e.g.. iHughes et al.lll994l ). However, as also noted 
by iNeininger et~aD Jl998h and IPetitpas fc Wilson! j200of )7 for optically thick radiation, it 



will be possible to directly detect emission only from the foreground surface of the imaged 
structure. In that case, the dust clouds will be seen as two emission peaks, while the torus 
(or bar) will manifest as an elongated, barlike emission of roughly the same optical depth or 
optically-thick brightness. 

Second, the two main peaks seen in maps of similar resolution as the 450 and 350 fim 
images in Figure s [J are not symmetric. In maps of better resolution and sensitivity (e.g., 
Shen fc Lol Il995l . see the CO contours in Figure [3]), the peak west of the nucleus has a 
morphology clearly different from the eastern peak and can be resolved into t wo or three 



struct ures. High-resolution maps obtained with the Very Large Array ( VLA) by IWills et al. 



(119991 ) showed that the western peak is associated with locations of supernova explosions 
of higher int ensity and e a rlier e volutionary stage than the eastern peak and confirmed the 
discovery by IWeiB et al.l (119991 ) of an expanding "superbubble" that is centered near the 
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location of M82's brightest supernova remnant, 41.9+58, and the submm western peak. 
The varying supernova intensities and ages across the disk of M82 is supported by high 
resolution HST imaging of stellar clusters that indicates that the re gions near the eastern 
and western submm peaks have different star-formation histories (e.g.. lde Grijsll200ll ). In this 
light, the submm peaks indicate concentrations of dust environments associated with different 
locations of varying supernova and star-formation activity, not the commonly assumed dust 
torus. 

Third, observations of line ratio gradients indicate that the average temperature across 
the lobe increases fro m the northeast to sout hwest, while the density increases in the op- 
posite direction (e.g.. iPetitpas &: Wilsonll2000l ). Further evidence of the higher temperature 



or, at least, column density in the southwest is seen in the lope-sided 850 and 750 /im lobe 
and the double-peak 450 and 350 /xm lobes in which the southwest parts are t he brighter. A 



torus that probably houses and is heated by an active galactic nucleus (e.g., iMuxlow et al. 



19941 ). would be expected to have a temperature that decreases from its inner to outer walls. 
One explanation for the higher temperature and lower density in the the southwest is linked 
with star-forr nation activity th at both heats and depletes the interstellar medium (ISM) at 
this location (IWills et al.lll999l ). or sim ply western and eastern reg i ons of two differen t star 



formation physical environments (e.g., lAchtermann fc Lacyl Il995l : iForster Schreiberl 12000 



de Grijsl l200ll : iLipscy fc PlavchanI |200J) . Evidence in this paper in terms of clearly asym- 
metric submm emission seen in almost all the presented intensity maps seems to disfavor the 
commonly assumed interpretation of a dusty torus in M82. 



5.2. Fluxes and Spectral Energy Distribution Analysis 

Section m above raised the possibility that CO emission may contaminate the continuum 
of M82 observed in SCUBA filters. For one, significantly different contributions in the 
SCUBA bands imply different corrections to the measured fluxes and would affect the SED 
analysis using those fluxes. If the differences are significant the SED computation and 
analysis should in theory only be determined after correcting or accounting for the CO 
contamination. It is currently not obvious if the contamination to the higher-frequency 
continuum will b e less or higher than the 4'^% estimated to the SCUBA-850 /im band by 



Seaquist fc ClarkI (120011 ): even though it is clear that it may also be important (see SectionH]). 
While relevant CO data need to be acquired in the future to make a quantitative calculation 
of the relative contributions in the SCUBA bands, for the practical determination of the 
SEDs of M82 in this paper, it is assumed that the CO contribution in the SCUBA bands is 
the same. 
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Table 2: SCUBA continuum fluxes for specific locations in M 82 (The listed measured flux densities 
include roughly 47% possible contribution from CO as discussed in the text.) 



Locations parallel to M82's 


Flux (Jy) 


Flux (Jy) 


Flux (Jy) 


Flux (Jy) 


disk position angle of 72° 


@350/im 


@450/im 


@750^m 


@850/im 


©Peak Flux 


18.6 ±5.6 


13.0 ±2.6 


1.8 ±0.4 


1.2 ±0.1 


30" X 15", about nucleus 


49.3 ±13.7 


28.8 ±5.6 


3.2 ±0.5 


2.3 ±0.2 


70" X 40", about nucleus 


63.4 ± 18.9 


35.9 ± 7.0 


7.7 ± 1.5 


3.8 ±0.4 


30" X 15", 15"N of nucleus 


15.9 ±4.8 


7.6 ±1.5 


1.0 ±0.2 


0.6 ±0.1 


30" X 15", 15"S of nucleus 


12.8 ±3.8 


8.1 ±1.6 


0.7 ±0.1 


0.6 ±0.1 


30" X 15", 30"N of nucleus 


6.2 ± 1.9 


2.8 ±0.5 


0.3 ±0.1 


0.2 ±0.02 


30" X 15", 30"S of nucleus 


3.1 ±0.9 


1.8 ±0.4 


0.2 ±0.1 


0.1 ±0.1 



Table [2] shows the submm fluxes measured at specific locations across M82 using 
SCUBA. The listed errors include calibtation uncertainties. The submm dominant emit- 
ting region in M82 is within 30" x 15" about the nucleus and is associated with the most 
intense star-formation in the galaxy. Assuming the primary source of the submm emission is 
dust re-radiation, the measured fluxes from this galaxy were fitted with the following thermal 

= nB,iT)[l - exp(-(^)^)], (1) 

A 

where Q is the solid angle for the emitting region, B^{T) the Planck f unction at temperatu re 
T, Ao the wavelength at which the optical depth is unity (Ao = 7.8 /im. lHughes et al.l ( jl994j )). 
and (3 the emissivity index of the grains. Due to the limited frequency sampling in the data, 
the temperature T and emissivity index (3 were the only parameters that were statistically 
determined in Equation [TJ 



Table 3: SCUBA derived dust emission properties for locations in M82 

Locations parallel to M82's T (3 Q 

disk position angle of 72° (K) (sr) 

30" X 15", about nucleus 36 ± 14 2.0 ±0.1 1.66e - 8 

30" X 15", 15"N of nucleus 31 ± 3 2.2 ± 0.1 1.66e - 8 

30" X 15", 15"S of nucleus 31 ± 15 2.3 ± 0.2 1.66e - 8 

70" X 40", about nucleus 27 ± 5 1.8 ±0.1 9.05e - 8 



Table [3] lists the derived dust T and (3, as well as the fl associated with the specified 
emitting regions in M82. For the listed rectangular locations, the Q is determined from 
the longest side of the rectangle. The average derived T is ~ 31 K, with a minimum and 
maximum of ~ 27K and ~ 36 K respectively. The highest T's are from regions of the highest 
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surface brightnesses, presumably corresponding to regions of intense star-formation, while 
the coldest T's are at regions farthest from the nucleus. The f5 values range from 1.8 to 
2.3 and are highest at regions farthest from the galactic nucleus. The higher (3 values are 
associated with larger and colder dust grains; therefore, the spatial SED analysis here points 
to colder grains more prevelant with increasing distance from the galactic center of M82, 
particularly along the minor axis of the galactic disc (see Section [573] below) . 

If the CO contamination in the SCUBA bands were not constant for M82, as assumed in 
this paper, an SED with SCUBA fluxes may have looked different than the result above and 
possibly required a different interpretation. For example, a CO contamination that increased 
with fre quency, as suggested from an analysis of lower-transition CO and CO(4-3) data of 
M82 by lGuesten et al.l (Il993l ) (see Section H]), would mean that SCUBA fluxes uncorrected 
for CO contamination lead to the SED indicating lower dust temperatures and/or higher 
emissivity indices than was actually the case. A proper correction for any CO contamination 
is therefore needed before a more definitive SCUBA SED analysis can be conducted for 
M82 and is worth attempting in the future, when CO data in all the SCUBA wavebands 
(especially the 350 and 450/im ones) are available. 



5.3. Implications of the SCUBA data and SEDs on the Outflow of Cold Dust 



In has been noted that the continuum morphology in the inner-halo of M82 has a 
general north-sout h asymmetry, that is consistent with the north-south asymmetric X-ray 
and Ha winds (e.g.. lWatson et al.lll984l . and see Section [3. 2. 31) and the associated UV, optical, 
molecular and irideed infrared to mm structures that have been reported in M82 (e.g., 
Seaquist fc ClarkI 120011 ). Therefore, a simple interpretation of the asymmetric morphology 
of the submm continuum in the halo of M 82 and is that it is a manifestation of an outflow 
of dust from the inner disk to the halo. 



Seaquist fc Odegardl (Il99ll ) presented one of the earlier, extensive evidence of the disk- 



to-halo outflows from spectral index distribution computations using radio continuum maps 
at several wavelengths between 0.33 and 4.9 GHz (90 and 6 cm). They found spectral indices 
between -0.3 and -0.6 in the inner-disk region, steepening to about -1.0 at a radius of about 
Ikpc along the minor axis, and concluded that these were from relativistic synchrotron- 
emitting electrons that were being scattered aga inst infrared photons emitted in the inner- 
disk region of M82. Recently, IWills et al.l (119991 ) used high resolution VLA continuum data 
between 1.4 and 5 GHz and computed spectral indices from —0.6 to —0.8 about 20" north 
of the disk, in localized nuclear sit es of the onflows tha t they call 'chimneys'. These values 
are consistent with the results by ISeaquist fc Odegardl (Il99ll ) in the same wavebands, and 
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Wills et al.l (119991 ) also interpreted them as indicating synchrotron emission from relativistic 
electrons entrained in the wind. 

A thermal component in the filaments has previously bee n suggested based on 'tentacles' 
observed in Nell maps of M82 (lAchtermann fc Lacyl Il995l ). which presumably are also a 
manifestation of the outfiow phenomenon. In comparative SED analysis of 30" x 15" regions 
centered in galactic nucleus and two others 15" north and south of it, the regions north 
and south had respectively cooler temperatures and higher emmisivity indices than the 
central region. This spatial SED analysis is consistent with the submm emission coming 
from a thermal source with a temperature decrease and emmisitivity index increase along 
the minor axes of the disk of M82. The change of the dust properties along the minor axes 
has a simil ar direction to the radio spectral index gradient shown by ISeaquist &: Odegard 
(Il99lh and Iwills et all imi\ and is con sistent with the north-south asymmetric winds 
(e.g., IShopbell fc Bland-Hawthornlll998l ) of M82 that has been shown to be co-spatial with 
the submm morphology (see Figure [5] and Section [3.2.3l) . In this light, the submm continuum 
maps (and SED changes along the minor axes) indicate an outfiow of dust grains that are 
ejected from the inner disk by, or entrained in, the starburst winds. 



One explanation for the outfiow being asymmetric was given by lShopbell &: Bland-Hawthorn 



(119981 ) ■ who suggested that if the star-forming disk is slightly shifted up from the galactic 
disk, this would imply that there is less covering material in the north and would make 
collimation difficult, resulting in an immediate blow-out of material in the north. Detections 
confirming that collimation is better to the south of M82 have been made of large-scale 
emiss ion extending to 1.5kp c, and more extended in the south, not only in opt ical line maps 
(e.g., Devine &: Ballvlll999l). but als o in C0(2-l) and CO(3-2) respectively by iThuma et al. 



(I2OOOI ) and lSeaquist fc ClarkI (I2OOII ). Another valuable result of the co-addition of SCUBA 
archive data in this paper is that the most sensitive submm maps that are displayed in 
Figure [T] show submm extended emission that is associated with the outfiows on scales that 
for the first time match the 1.5 kpc CO detections noted above. 



Recently reported radio 'chimneys' (e.g.. IWills et al.lll999l ). that are about 20" north 



of the disk and hypothesised to signify local blow-outs of material by supernova-driven 
winds, are not obvious in the SCUBA maps. However, prominent SiO features associ- 
ated with the localized radio outfiows have now been detected in mm-heterodyne obser- 
vations obtained with the Institut de Radio Astronomic Millimetrique interferometer (e.g., 
Garcia-Burillo et al.ll200ll ). These authors explain the SiO detections in a framework of 
shocked chemistry at the sites of the gas ejections from the starburst disk. For the moment 
it appears that the shocked gas has proved to be a better probe of the localized onflows than 
the direct observations of dust emission in the mm to submm continuum. 
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Localized sites of the outlows or 'emission spurs', althoug h long sought-after and some- 
times reported in mm to submm continuum and CO maps (e.g.. lHughes et al.lll994j : IShen &: Lo 
19951 : iLeeuw et al.lll999l : lAlton et al.lll999l ). have not been reliably reproduced in the different 
mm to submm observations. All the continuum maps in this paper also have some low-level 
'spurs', but almost none are reproduced at exactly the same locations and extend to the same 
degrees between any two different observations. This would suggest that the spurs might 
be artefacts in the maps. However, if the mm and submm spurs are emission from dust 
outflows (and filaments, clouds or lanes) that are of different optical depths or compositions 
and re-radiate low-level emission in relatively narrow wavebands, the submm spurs may not 
be reproduced in maps at certain wavelengths and sensitivities. 

It was noted in Section 13.21 that dust in M 82 is most probably within components of 
physically varying locations and origin or simply at different spatial-depths or stratifica- 
tions. The r e cent high-resolution, interferometry observations of M82's molecular gas by 
Walter et al.l (120021) discovered CO (1-0) 'streamers' and decoupled previously ob served CO 



outflo ws (e.g.. lNakai et al.lll987l : ISeaquist fc Clarkll200ll ) from the streamers (e.g.. lYun et al. 



19931 ). clarifying the spatial distribution and origin of the molecular gas in this galaxy. Sim- 
ilar, future observations with recently commissioned sensitive and high-resolution submm 
imaging instruments such as Spitzer and the SMA should provide tighter constraints on the 
spatial and optical depth properties of dust, test the reality of the reported submm spurs, 
and possibly detect dust re-radiation streamers, clarifying the implications or associations 
of all these features to dust outflows and recychng in M 82 (c.f. Section | 



6. Summary of Results and Future Work on M 82 

SCUBA 350, 450, 750 and 850 fim imaging observations have been presented of the 
dust-laden, star-forming inner disk and large-scale, low-level emission that is associated with 
the outflows in the halo of M82. The displayed maps include co-added data that were 
mined from the SCUBA Data Archive, resulting in the deepest submm continuum maps of 
M82. The 850 //m morphology has a single emission peak that is centered about 9" west 
of the galactic nucleus, while the 450 and 350 /im maps have two emission peaks centered 
about 10" and 6" respectively east and west of the nucleus along the galactic disk, similar to 
previous continuum observations at mm, submm and mid-infrared wavelengths, as well as 
to CO line transitions and Hq, observations of comparable resolution (see Section [3]). In the 
750 /im image, the single emission peak seen in the 850/um image (see Figure [1] is predictably 
beginning to be resolved-out into the double peaks seen in the 450 and 350/im images. Low- 
level emission is detected out to 1.5 kpc for the flrst time in the 850 /xm continuum of this 
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galaxy, i.e. at least ~ 160 pc radius farther-out than other recent studies. 



The deep maps were used in a detailed morphological study of the disk and large- 
scale detections, including a comparative analysis of submm to optical morphology (see 
Section [3.21) . The overall, extended submm morphology of M82 generally resembles the 
optical picture in that the disk emission has an apparently elliptical shape whose major axis 
is clearly aligned with that of the galactic disk at position angle ~ 72°. However, the submm 
morphology is much smoother than the optical picture, and some prominent dust cloud and 
filamentary lanes that are seen in the optical are not obvious in the submm continuum. One 
simple explanation for the submm versus optical correspondence (or lack of it) is the different 
resolutions and sensitivies of the presented observations. If the optical features emit submm 
emission, it is possible the emission is at a level lower than the current mapped SCUBA 
sensitivities or smeared and thus erased by the relatively larger SCUBA beam. This can 
be varified by future higher resolution and more sensitive submm observations that should 
become possible with ALMA. 

A comparative analysis was further conducted of submm to high-resolution CO (1-0) 
morphology (see Section 13. ip . Some resolved peaks in the CO maps could be associated 
with unresolved features in the submm maps. However, there are differences that show the 
CO and dust emission do not trace each other in a very simple way and suggest that although 
the dust and CO generally appear mixed in the central star-forming region, in fact there 
are differences in their spatial distributions. The different concentrations are most probably 
locations of varying ga s-to-dust densi t ies or star-formation env ironments, that have been 
reported in M82 (e.g., Ide Grijd l200ll : iPetitpas fc WilsonI 120001 ). This can be varified by 
future higher resolution observations of dust in M 82 that should be possible with the SMA 
and later ALMA. 

The SCUBA maps were also used in a computation of the first submm spatial SED 
analysis of locations within and outside the central star-forming region of M82 (see Sec- 
tion [5]2]) and in the discussion of the origin and structure of submm maps (see Section [5]). 
In particular, (a) the commonly assumed interpretation that the double emission peaks that 
were seen in the mm to infrared continuum are due to emission from the edges of an in- 
clined, dusty-molecular torus was challenged (see Section [S.ip . (b) an analytical review of 
CO results was undertaken to assess if CO emission might significantly contaminate the 
continuum observed in SCUBA filters (see Section H]), and (c) a morphological comparison 
was conducted to ch eck whether the lo caliz ed outflows that were repor ted in radio and SiO 
maps respectively by lWills et al.l (119991 ) and lGarcia-Burillo et al.l (120011 ) could be seen in the 
SCUBA maps (see Section [5T3|) . 



Evidence in this paper in terms of clearly asymmetric inner-disk submm emission seen 
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in the presented intensity seems to disfavor the commonly assumed interpretation of a dusty 
torus in M82. Arguments were presented to explane the inner-disk submm maps of M82 
in the context of emission from a rather complex distribution of dust concentrations that 
are in regions of different star-formation envi ronments, as has been reported from various 
studies using data a t other wavelengths (e.g., lAchtermann fc Lacylll995l : iForster Schreiber 
2000l : lde GriSboOlh . 



It is not obvious if the CO contribution to the higher frequencies of S CUBA will be less 
than or higher than the 47% estimated to the SCUBA-850 fim band by ISeaquist fc Clark 
(I2OOII ). However, it is clear that the CO contamination to higher-frequency continuum may 
also be significant and warrants detailed investigation. Future work on this galaxy will 
attempt to acquire the data of the CO lines in the 450 and 350 /im bands and make a 
quantitative comparison of these data in order to determine the possible contributions of 
CO to the high-frequency-SCUBA data. 

The overall submm low-level morphology has a general north-south asymmetry that is 
similar to the H„ winds and CO a nd X-ray outflows that have been detected in M82 (e.g., 
Shopbell fc Bland-HawthornI Il998l ). The submm spatial SED analysis also shows thermal 
properties that ch ange along the minor axis of the galaxy disk, s imilar to the radio spectral 
index gradient by ISeaquist &: OdegardI (Il99ll ) and IWills et al.l (Il999l ) and consistent with 
the north-south asymmetric, large-scale X-r ay and H„ winds. Therefore, the cu rrent re- 



sults support the simple interpretation (e.g., iLeeuw et al.lll999l : lAlton et al.lll999l ) that the 



asymmetric morphology in the submm maps is a manifestation of corresponding outflows of 
dust grains from the galactic disk into the halo. As noted above, this work has presented 
low-level 850 /im continuum emission out to 1.5 kpc for the first time in this galaxy, i.e. at 
least ~ 160 pc radius farther-out than other recent studies. 
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Fig. 1. — Deep maps of the 850 (top left - a), 750 (top right - b), 450 (bottom left - c), 
and 350 /xm (bottom right - d) continuum emission centered at the near-infrared nucleus of 
M82. The 850, 750, 450, and 350 yum contours on the respective maps are [16, 32, 48, 64, 
96, 128, 200, 300, 400, 600, 800, 1000, 1175, and 1350], [200, 400, 600, 800, 1000, 1300, 
1600, and 1900], [450, 900, 1450, 2000, 3000, 4000, 5000, 6000, and 6800], and [1300, 1950, 
3250, 5200, 7000, 8500, and 10000] mJy/beam; the respectively rms uncertainties are ~ 8, 
100, 225, and 650mJy/beam, and the resolutions are ~ 14".5, 11".4, 8".5, and 6".7. The 
white circle in each map indicates the approximate size of the SCUBA beam at the plotted 
wavelength. The keys are grayscale-coded intensities in Jy/beam, and the X and Y axes are 
J2000 coordinates. 
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Fig. 2. — Deep maps of the 850 fim continuum emission (left - a) and integrated C0(2-l) line 
intensity ( right - b) of M 82 res pectively as presented in Figure [T] of this paper and Figure 5 
of work by iThuma et al.l (120001 ) . The respective beams for the two observations are plotted 
and the rest of the keys, contour lines, axes are as presented in the indicated figures. A scale 
bar of ~ 1.5 kpc is shown on the maps. 
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Fig. 3. — The central 76" x 56" region of M82 at 850 (left panels - a & c) and 450 /im 
(right panels - b & d) overlaid with the i ntegrated-intensity (black contours) of the CO (1-0) 
interferometry data by lShen fc Lol (119951 ). The top panels are plotted at full resolution of the 
CO contours and the bottom panels are with the CO contours smoothed to beam sizes similar 
to the SCUBA-850 and -450 /xm beams (~ 14". 5 and 8". 5, respectively). The white 850 and 
450 /im contours on the respective maps are [200 to 1175] and [900 to 6000] mJy/beam, as 
displayed in Figure [H The red 850 and 450 yum contours are at 1350 and 6800mJy/beam 
and are plotted to highl ight the position of the respective western submm peak. The cross 
marks the 2.2 //m peak (jPietz et al.lll986l ). The keys are color-coded intensities in Jy/beam 
and the X and Y axes are J2000 coordinates. 
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Fig. 4. — The left (a) and right (b) p anels respectively sho w the B-band maps of M82 
obtained by the Hubble Legacy Team (IMutchler et al.l 120071 ). respectively overlayed with 
the 450- um continuum-era ission contours shown in Figures [1] and CO (1-0) interferometry 
data by IShen fc Lol (119951 ) . For best contrast, the optical intensities are inversely plotted 
and therefore the light patches are extinction feat ures. The centeral cross and open squares 
respectively mark p ositions of the 2.2 nm peak by iDietz et al.l (119861 ) and mid- infrared, star- 
forming clusters by lLipscy fc Plavcharu (120041 ). The X and Y axes are J2000 coordinates. 
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Fig. 5. — The left (a) and right (b) panels respectively show t he low- and high-brig htness 
features of the Hq, maps obtained by the Hubble Legacy Team (IMutchler et al.ll2007l ). over- 
layed with SCUBA 850 and 450-/im continuum-emission contours as shown in Figures [H 
The overlays were made by aligning the sky co-ordinates of the Hubble maps with those of 
the SCUBA ones. The respective resolutions in the Hubble, SCUBA 850-//m, and 450-^m 
maps are < 1", ~ 14". 5, and ~ 8". 5. The overlay in the left panel demonstrates the spa- 
tial co-incidence between the large-scale low-brightness Hq, emission and 850-fim continuum, 
while the right panel depicts the origin of the large-scale H^ emission in the intensive star- 
formation inner region of M82's disk about the submm emission peaks. For best contrast, 
the optical intensities are inversely plotted and the light patches (e.g. across the center of 
the image in the right panel) are foreground extinction features. The X and Y axes are J2000 
coordinates. 
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